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We consider a model in which dark matter is stable as it is charged under a Z2 symmetry that is 
residual after an A4, flavour symmetry is broken. We consider the possibility to generate the quark 
masses by charging the quarks appropriately under A4. We find that it is possible to generate the 
CKM mixing matrix by an interplay of renormalisable and dimension-six operators. In this set-up, 
we predict the third neutrino mixing angle to be large and the dark matter relic density to be in 
the correct range. Low energy observables - in particular meson-antimeson oscillations - are hard 
to facilitate. We find that only in a situation where there is a strong cancellation between the 
Standard Model contribution and the contribution of the new Higgs fields, B meson oscillations are 
under control. 

PACS numbers: 

I. INTRODUCTION 

Nowadays, there is strong observational evidence of the existence of dark matter (DM) [TJ [3J. Many experiments 
are currently looking for direct or indirect observation of a dark matter candidate [3H5]. Among all the possible 
DM candidates stable cold dark matter (CDM) ones have been discussed in many Standard Model (SM) extensions. 
Typically its stability can be secured by introducing a parity, under which the CDM candidate is odd. This parity is 
often introduced ad hoc with the sole purpose of making the DM candidate stable or it is related to a parity added 
to the theory by hand. 

Recently, a model [SHE] (to be referred to as discrete dark matter or DDM) was proposed that relates this DM 
parity to the residual symmetry of a spontaneously broken flavour symmetry. Flavour symmetries became popular 
after the discovery of particular patterns in neutrino mixing, since they can reproduce the observed structures from 
symmetry principles. 

In DDM, only the lepton sector was studied. In this work, we will consider a simple way to add quarks to the 
model. We will show that this results in a diagonal CKM quark mixing matrix at the renormalisable level, but that 
non-renormalisable operators can generate correction to this in order to reproduce the observed mixing patterns. It 
is not the first time that a discrete lepton non-Abelian flavour symmetry is extended to the quark sector [5]-[2"5]. 
However, contrary to the majority of the cases in our scenario we do not have flavons, that is heavy flavour scalar 
SM singlets; instead SM scalar doublets transform non trivially under the flavour symmetry. A consequence is the 
appearance of multiple new Higgs fields. Recently a similar set-up has been proposed in |29j . 

We will see that in our scenario we can predict the third neutrino mixing angle to be in the near-future experimental 
sensitivity (and indeed in the T2K range [30 ) and that we have a strong enhancement of meson-antimeson oscillations. 
These are so strong that they rule out almost all parameter space for the model. Only if we allow a strong cancellation 
between the well-known contribution of the Standard Model and new contributions AMg s , AMg d and AMx can be 
reconciled with the data. Nevertheless we see only a limited effect on the calculation of the relic density of the DM 
candidate. 

II. THE MODEL 

As in DDM, we assign matter fields to irreducible representations of A4, the group of even permutations of four 
objects, isomorphic to the symmetry group of the tetrahedron. The group A4 is often used in flavour model building 
as it naturally allows the neutrino mixing to be of the tribimaximal type |31j . that fits the observed data very well. 
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The properties of the group A4 are summarized in the appendix. 

The representations of the leptons and the Higgs fields are as in DDM. In particular, there is an Sf7(2)-doublet 
Higgs H in the trivial representation of A4 and we assume three extra copies of the Higgs 77 = (771, 772, 773) transforming 
as a triplet. There are also four neutrinos: three of them transform as an ^-triplet Nt = (Ni, N 2 , N 3 ) and a fourth 
transforms as a singlet N4. We assume that all quarks transform in the same way as the charged leptons. Both 
lefthanded and righthanded fields transform as one-dimensional representations of A4, with different representations 
over the generations: the first generation is taken to transform as 1, the second generation as 1' and the third 
generation as 1". They thus transform non-trivial under the Z 3 subgroup of A±, but are uncharged under the Z 2 
subgroup. All matter and Higgs assignments of our model are summarized in table]!} 
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TABLE I: Summary of relevant model quantum numbers, q = u,d. 



The resulting Yukawa Lagrangian for the leptons is unchanged 



A = y e LJ Re H + y^L^r^H + y T L T l c Rr H+ 

y^MNrfj^ + y^L^Nr^r, + y v 3 L T (N T fj) v + (1) 
y%L e N A H + MxN T N T + M 2 N A N A + h.c, 

where we have defined 77 as iu 2 rf . 
For the quarks, the Lagrangian reads 

£ q = VuQiHuxr + y c Q 2 Hu 2R + y t Q 3 Hu 3R + 

ydQiHdiR + y s Q 2 Hd 2 R + ybQ^dzn + h.c, 

where H stands for H, defined analogously to 77. 

At the renormalisable level, both up- and downquark matrices are diagonal with all masses given by m; = yiVji / v2j 
with Vh/\2 the vacuum expectation value of the A4-singlet Higgs field. As we do not aim to explain the hierarchy of 
the quark masses, with m u 10~ 5 m t , we assume an hierarchy in the Yukawa couplings, with most Yukawa coupling 
being small to very small. For instance a Froggatt-Nielsen symmetry |32j could make this more natural. 

As discussed in more detail in [B], electroweak symmetry is broken by the vacuum configuration 

(if ) = v H /V2, (77?) = (v v /V2, 0, 0) . (3) 

We write the ratio between the vev of H and 77 as tan/3 and obviously, their squares sum to (246 GeV) 2 . The vev of 
77 breaks the A4 group into its subgroup Z 2l generated by S, that is diagonal in the three-dimensional representation: 
S = Diag(l, —1, —1). The Z 2 symmetry thus acts on the A4 triplet fields in the following way: 

z N 2 -> —N 2 , h 2 -> -h 2 , A 2 -> -A 2 , , 4) 
2 ■ N 3 -> ~N 3 , ha -> -h 3 , A 3 -> -A 3 , 

where ^2,3 are the components of the triplet Nt and ft.2,3 and ^2,3 are respectively the CP-odd and CP-even compo- 
nents of the Higgs doublet 772.3. 

The residual Z 2 symmetry is responsible for the stability of the lightest combination of h 2l h 3 , A 2 and A 3 which 
is the dark matter candidate. Indeed the Z 2 -odd candidate only couples to heavy right-handed neutrinos and not to 
the SM charged fermions, that are Z 2 -even. Such a scalar dark matter candidate is potentially detectable in nuclear 
recoil experiments [H |S] . 

We refer to the four i?2-even components of the i?2-even Higgs fields H and 771 as Hq, H[, A' q and A\ in accordance 
with [8]. They give rise to two scalars H and Hq, one pseudoscalar and the neutral Goldstone boson of electroweak 
symmetry breaking. 
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As mentioned above, the charged leptons and quarks transform non-trivially under the subgroup of A4 generated 
by T. The vev configuration ^ clearly breaks this Z 3 . Still, at the tree level, quark and charged lepton masses preserve 
the Z$ invariance thanks to the scalar charge assignments. As we will see in the next section Z3 breaking effects appear 
at next to leading order (NLO) level giving rise to the quark mixing matrix. 



III. QUARK MIXING 

Quark masses at the tree level are given by equation ([2]). This gives rise to diagonal quark mass matrices and the 
CKM matrix Vckm = (V^yV^ is simply the identity matrix. 

The gauge and flavour charge assignments in the Higgs sector allow the construction dimension six operators for 
the down-type quark masses that contain H, 77 and their conjugates. There are three ways to contract the SU(2) 
indices, represented by brackets in the equation below 



E 



|f (Q.i? K fat,,) + f JL (Qtrfdj fatF) + f ~§ (QiV)dj (frr,) ■ (5) 



The contraction of A4 indices between the two j/t) triplets (see equation (A3|) is such that it generates the right 
type of singlet (f, 1' or 1") to match the charges for Q i and dj- It is important to note that this is possible for any 
combination of i and j due to the product rules of A4. A is the cut-off scale, up to which we accept the theory to 
be valid and the / couplings are dimensionless. Analogous dimension-6 operators can obviously be constructed for 
up-type quarks and charged leptons. 

The mass term Lagrangian ^ and the effective couplings ([5| generate the effective mass matrix for down-type 
quarks 

m d \ v v 2 ( h dd h ds h db 
-U;= j m s +^ h sd h ss h sb I + Od/A 1 ). ((,) 
m b ) \h bd h bs h bb 

where fty = (/ y + /(, + f 1^/2^2. Analo gous expressions again hold for up-type quarks and charged leptons. 

Now the crucial question is how large the cut-off scale A is. In principle, this is a scale we are free to set. Only 
using 'naturalness' and 'finetuning' arguments, we can find a range for it. We will give two arguments, both pointing 
to a scale of I to 10 TeV. 

In the first argument, we demand that there should not be more than 10 to 100% corrections to the Higgs from 
one-loop corrections to the Higgs propagator with the fermions and the (new) scalars of the theory. These corrections 
are typically of the order A 2 /(47r) 2 and requiring them to be not too large with respect to v^ w — (246GeV) 2 indeed 
gives A < (1 to 10) TeV. 

Interestingly, we find the same scale from an argument where we require the dimensionless parameters h, in particular 
h ds , to be of order 1. The off-diagonal terms in equation ^ are responsible for generating the quark mixing as 
parameterized by the CKM matrix. As we do not have information about the size of the dimensionless parameters h, 
we assume them to be of order 1, which can be seen as the most natural assumption for dimensionless parameters. 

Under this assumption, the absolute values of the corrections to the leading order elements of the mass matrix are 
of the same order for the up-type quark matrix and the down-type quark matrix. However, due to the much larger 
elements of the leading order up-type quark mass matrix, the effects on quark mixing are dominated by the down-type 
quark contributions. This allows us to estimate the order of magnitude of the cut-off scale. 

Now the (1 2) element of equation Q should be of order \cra s in order to reproduce the Cabibbo angle. 

jY 2 " = *cm s - (7) 



This gives 



Xcm s (tan 2 /?)(! + ^) 3/2 A c ™ s 
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depending on the exact values of hds and tan/?, which we have taken between 0.1 and 1 and between 0.1 and 10 
respectively. Due to the large bottom mass, the effect of hdb and h s b of the same size as hds on #13 and #23 is relatively 
minor and these angles are thus naturally smaller than the Cabibbo angle and indeed in a large part of parameter 
space, we can fit them to their measured values. 

The analogue of the dimension 6 operator ^ affects the lepton mixing. In [5], it was shown that at leading order, 
the lepton mixing matrix has zero #i3-angle and large 9±2 and #23 mixing angles, although these do not necessarily 
fit in a mixing pattern such as tribimaximal or bimaximal mixing. The fact that the down-type quark and charged 
lepton mass matrices are alike (at least at leading order) suggests that the matrices that diagonalize them, V* and 
V£ are also similar. We thus expect a large angle (of the order of the Cabibbo angle) in the (1 2) sector of V£ and in 
the lepton mixing matrix (V£)^ V£, this affects all thr ee angles. In particular, we expect a Cabibbo-sized correction 
to the #i3-angle. This is in the range of values hinted at by the recent T2K results [30] and in recent global fits [55] , 
In any case, we predict that non-zero #13 should clearly be measured by next generation experiments, such as Daya 
Bay [35 and Double Chooz [35]. 

We end this section with a comment on the scale of the neutrino seesaw. In the DDM model, neutrino masses are 
assumed to originate from the type-I seesaw 

m v = -m D3x4 M R ^ml 3 

This has two non-zero eigenvalues, that are of the order y\ v^/Mk- There can be a low energy seesaw, where 
we identify the mass scale of the righthanded neutrinos to A if the neutrino Yukawa couplings are not too small 
y v = 10 - ' 4 ^ 5 '. Lastly, next-to-leading order effects lift the mass of the lightest neutrino away from zero. Still, it is 
suppressed with respect to the other neutrino masses by a factor v^/A. 

IV. FCNC 

The inclusion of the operators eq. ([5| gives rise to tree level FCNC processes mediated by the ^-even scalar (Hq 
and H) and pseudoscalar (Aq). In the mass eigenstate basis, the trilinear couplings read 

l/(2V2)F Li f Rj Hv^v h ([U HHo } 11 /v h + 2[U HHa ] 21 /v v )(f lJ + ft + ft) , 
l/(2V2)F L J Rj H v 2 n v h ([U H H Q }i2/v h + 2[U HHa } 22 /v n )(f lJ + ft + ft) , 

±l/(2V2)F Li f Rj A [U GAo } 12 (f ij + ft - ft)/v h + 2[U GAo ] 22 ft/v v , (10) 

where in the last equation + is for up-quarks and — for down-quarks-and charged leptons as well. Uhh an d Uga 

are the matrices that relate the scalar and pseudoscalar mass and interaction eigenstates; in particular Uga is a 

rotation matrix over the angle j3 [8]. 

In this section, we will focus on meson-antimeson oscillations, as these are among the most constraining tests for 

new physics. In particular, we find that they are more constraining than the often-discussed K-meson decays. In 

our model, the Standard Model box diagrams are accompanied by new tree diagrams; see figure [T| In models with 

— 

multiple Higgs bosons, the new physics (NP) contribution to the mass splitting connected to F° — F oscillations is 
given by [HUH] 

at 

Here, Mp is the mass of the meson, fp is its decay constant and Bp are recalibration constants of order 1, related 
to vacuum insertion formalism. The masses, m r and m s are those of the quarks of which the meson is build, i.e. 
rs = bd, bs, ds stands for Bd, B s and K° respectively. Lastly, I" s and J" s are effective fermion-fermion-scalar and 
fermion-fermion-pseudoscalar couplings, as given in [38] . 
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FIG. 1: Feynman diagrams for B mesons oscillations in the SM (left) and in our model (right); similar diagrams can be drawn 
for B a , and K mesons. 



Indeed we find that the meson-antimeson oscillations severely reduce the parameter space of our model. As men- 
tioned in the previous section, the CKM matrix can dominantly originate from corrections to the up-type quark mass 
matrices or the down-type quark mass matrices and we mentioned that dominance of the latter is more natural. 
Indeed we find that if corrections of the former type dominate or even if there is no dominance of one of the two, 
AMrj of D meson oscillations is much larger than the experimental value [37] and this scenario should be excluded 
as shown in figure [2] 
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FIG. 2: D meson oscillations in a scenario where the CKM matrix is generated dominantly by corrections to the up-type quark 
mass matrix. 



In case the CKM matrix mostly originates from corrections to the down-type quark mass matrix, we find that the 
bounds from meson mixing are rather strong. In figure |3j we show the contribution of the new diagrams to these 
as well as to AMs d as function of the lightest Higgs mass. We see that in almost all of parameter space the points 
are near or even above the short-dashed line, which indicates the current experimental value 39J that is rather well 
described by the Standard Model box diagrams [101 SI] • Naively, this is interpreted as an exclusion of the model, 
which is true in most of parameter space, but not in points where the Standard Model and new physics contributions 
partially cancel. To see this, we write the mixing amplitude for Bd mixing as [32] 

Mnp = a - 1+ ^j [mu^ (12) 

and analogously for B s . We have verified that our expression for the NP contribution carries enough phases to 
generate a flat distribution for ad and a s . We check that the points in (h d = 0.41, ad = 100°) and {h 2 s = 1.6, a s = 90°) 
are allowed by the data [42] [43] and give a nett contribution of NP with respect to the full result of respectively 0.65 
and 2.67 times the observed value. For Bd mesons, NP effects are forced to be less than the observed value, while 
for B s mesons, it can be slightly more. These values (and a corresponding estimate in case of the kaons) correspond 



() 



to the dot-dashed lines in [3] We see that a small, but significant number of points is allowed by the data under the 
assumption of partial cancellation. 
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FIG. 3: B and K mesons oscillations. The horizontal lines are as explained in the text. 

The naturalness of this cancellation, that requires hnctuning between the phase and amplitude of the new and 
the Standard Model contribution, can be a matter of debate. The need for cancellation diminishes for larger Higgs 
masses, although there are still no points below the lower line for B s and K. Indeed, requiring that the new diagrams 
contribute less than the experimental bound, as is customarily done, the model would be excluded. 

On the other hand allowing a strong negative interference between the SM and the DDM contributions does not 
further constrain the scalar spectrum with respect to the analysis done in [S] . Indeed fig. [4] shows that there is no 
correlation between the bound imposed and the mass of the lightest iSVeven scalar state, even if the number of points 
allowed significantly reduces with respect to those in [5]. 

V. DARK MATTER RELIC DENSITY AND DIRECT DETECTION 

The dark matter candidate of the DDM model can annihilate into up or down quarks via SM Higgs exchange in 
the s-channel. The operators ^ give rise to a new contribution for this decay in the case of decay into down quarks. 
Decay to up quarks is not affected as we have seen that up-type quarks are forced to be almost diagonal. In [5] it was 
shown that the model -without NLO terms - gives rise to an available DM candidate. The relic dark matter density 
constrains the parameter space of the model. In particular it was found that dark matter mass Mjjm is in the range 1 
- 100 GeV. The region with M DM < 40 GeV and mass of the Standard Model Higgs M H > 400 Gev is excluded, while 
for Mom > 50 GeV the Higgs mass can go up to about 500 GeV as co-annihilation can be possible. For Mdm lighter 
than about 80 GeV the dark matter annihilate (coannihilate) into fermions trough exchange of scalar (pseudoscalar 
and gauge boson) in the s-channel. For masses heavier than about 80 GeV the main channels of annihilation are with 
W bosons in the final state. 
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FIG. 4: The mass of the lightest Zi-even state versus the mass of the DM candidate without (red) and with (blue) meson 
oscillation constraints. 



Below we consider the effect of the inclusion of the dimension-six terms of eq. Q to the relic density. The oper- 
ators ([5]) give an effective quartic coupling of the dark matter with quarks. We therefore study the effects of such a 
operator only for dark matter mass below 80 GeV as for heavier dark matter mass the main annihilation is into W 
bosons and not in quarks. 

We recall that the DM candidate is one of the four neutral states (2 scalar and 2 pseudoscalar) arising by the mixing 
of the neutral components of 772 and 773 that in [5] were indicated as H' m + iA' m with m = 2,3. From eq. ([5| we derive 
the H' m , A' m coupling to fermions. For the couple H' m H' m with m — 2.3 the coupling is proportional to 

v H (fij + /(,- + W/A 2 = vnh^/A 2 , (13) 

while for A' m A' m it is proportional to 

M/^ + iW^M 2 - ( 14 ) 

If the DM candidate turns out to be the CP-even ^-odd lightest state we may estimate the contribution of the 
new operator to the total c xx _ +d . ( j.. Similar conclusions would be obtained considering the case in which the DM 
candidate is the CP-odd Z2-odd state. In the previous section we have seen that 

, vi 

h^VH—^nisXc, (15) 

for any couple (ij) 7^ 11, since for the first family the is required to be smaller to fit down mass. There we need 

v 2 

h d u v H -^~m d . (16) 
We define the A'j^ as the parameter of the four points interaction xx^idj- We can estimate its size as 

\ d - . h d y H _ , m s^c ( 7 , 



We compare the effect of this operator on the the av re i of the process XX ~> d%di to the effect of s-channel exchange 
of the SM-like lightest CP-even scalar. The SM-like Higgs couples to fermions proportionally to y q ~ m q /vH and it 
turns out that the new contribution is negligible if 



m s X c ^ m q A H 

2 < T ■ \ L °> 



In eq. ( 18 ) my l is the mass of the lightest CP-even neutral scalar and Ah ~ vw is the dimensional coupling that 
controls the interaction of the dark matter with the Higgs doublet H\X- Since v v ~ Vjj ~ ~ Ajj ~ v^y the new 
contribution is naturally subleading for the second and third generation. For the first generation, on the contrary the 
new contribution to XX ~~ ^ dd is of the same order as the old one. We can (conservatively) estimate that this channel 
is negligible only for for values Mdm > 1 GeV. 

NLO terms also give rise to non-diagonal scattering XX ~* didj, which was not included in the analysis done in 
[5]. In particular the scattering XX ~> db can be non- negligible for DM masses around a GeV. For this reason we 
conclude that the previous analysis is not affected in the range Mjj m > 5 GeV. We take this lower bound as a further 
constraint on the scalar sector parameter space, postponing a complete new analysis to the future 

Direct detection is not affected at all by the NLO terms: the quark flavour diagonal scattering contribution are 
subdominant with respect to those mediated by the scalar H while the off-diagonal one could only give rise to processes 
that are not kinematically allowed such as x + N ~~ * X + ■A/' + tt + + e~ + v, with AT a nucleus in the detector bulk. 



VI. CONCLUSIONS 



Neutrino mixing might be explained by a discrete non-Abelian flavour symmetry. If this symmetry is dynamically 
broken only in one direction, a residual symmetry survives. It is interesting that this residual symmetry may be 
responsible for the existence of a stable dark matter candidate. It has already been shown that this set up can 
describe the physics in the lepton sector and the dark matter abundance rather well. 

It is a natural question to see if the quark sector can also be described in such scenarios. In this paper we investigated 
this possibility in a particular model [SHE] • We found that if we add quarks in the same non-trivial representations of 
the flavour symmetry as charged leptons, the interplay between renormalizable operators and dimension-six operators 
can generate a realistic CKM matrix in a large portion of parameter space. It is possible to let the CKM matrix 
dominantly originate from the up-type quark sector as well as the down-type quark sector, although the latter case is 
more natural. 

The new dimension six operators should also be present in the lepton sector. As a consequence we have shown that 
the predictions of the original model, #13 = and = 0, are shifted away from zero, in the former case into the 
near-future observable region. At the same time for what concerns the quark sector the new operators can lead to new 
channels for flavour changing neutral currents and we have analysed their effects on meson-antimeson oscillations. We 
found that D meson oscillations rule out the scenario where the CKM matrix originates mostly from the up sector 
and that K and B Sj( i meson oscillations are much enhanced as well, with some of their amplitudes at least as large 
as corresponding the Standard Model amplitudes, bringing the model close to being excluded, although a certain 
number of points in parameter space can still be reconciled with current observations. If this is the case the scalar 
sector behaves exactly as described in [8]. 

Although these observations do not rule out the extension of the discrete dark matter model with quarks, severe 
finetuning is needed and honesty forces us to say that this 'natural extension' is less natural than we hoped. 
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Appendix A: The group A4 

All 24 elements of A4 are generated from two elements S and T with S 2 = T 3 = (ST) 3 = X. A4 has four irreducible 
representations, three singlets 1, 1' and 1" and one triplet 3. 

We can choose a basis in which S and T can be represented as (1, 1), (1,0;), (1,cj 2 ), with u> = e 27 ™/ 3 for the three 
one-dimensional representations and 
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for the three-dimensional representation. S and T on themselves generate the two maximal subgroups of A4, the 
Abelian Z 2 and Z 3 . The A4 multiplication rules are given by 

1 x r = r for all representations r, 

1' X 1' = 1", 1" X 1" = 1, 1' X 1" = 1, 

(A2) 

1' x 3 = 3, 1" x 3 = 3, 
3 x 3 = 1 + 1' + 1" + 3 + 3. 

Representing the two triplets as a = (ax, a 2 , a 3 ) and b = (61, 62, ^3), the elements of the last product are 

(ab)i = a ± bi + a 2 b 2 + a 3 b 3 ; 

(ab)y = aibi + uja 2 b 2 + ui 2 a 3 b 3 ; 

(ab)i» = ai&i + uj 2 a 2 b 2 + u>a 3 b 3 ; (A3) 

(o&) 3l = (a 2 6 3 , 0361,01^2) ; 

(ab) 32 = {a 3 b 2 ,aib 3l a 2 bi) . 
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